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CZ3 , Abstract -We analyze the dynamics of desorption of a polymer molecule which is pulled at one 

of its ends with force /, trying to desorb it. We assume a monomer to desorb when the pulling 
force on it exceeds a critical value f c . We formulate an equation for the average position of the 
n th monomer, which takes into account excluded volume interaction through the blob-picture of a 
polymer under external constraints. The approach leads to a diffusion equation with a p-Laplacian 
for the propagation of the stretching along the chain. This has to be solved subject to a moving 
boundary condition. Interestingly, within this approach, the problem can be solved exactly in the 
' trumpet, stem-flower and stem regimes. In the trumpet regime, we get r = ron^ where nj is the 

number of monomers that have desorbed at the time r. to is known only numerically, but for / close 
to f c , it is found to be r ~ / c /(/ 2/3 - fc /3 )- If one used simple Rouse dynamics, this result changes 
to r ~ fc n d/ (/ — fc)- I n the other regimes too, one can find exact solution, and interestingly, in 
all regimes r ~ n\. 
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Introduction. — Single molecule experiments, in which one exerts an external force 
on a long chain molecule have become ubiquitious pQ. These experiments provide very 
interesting information on the structure and dynamics of long chain molecules. Mechanical 
force can be used to unzip the double stranded DNA, emulating the action of enzymes 
that mechanically force the two strands apart during the process of replication [2J. Force 
has also been used to desorb a molecule adsorbed on a surface 00]. Motivated by the 
experiments, the statistical mechanics of an adsorbed polymer molecules, at one end of 
which a force is acting trying to desorb it has been extensively studied [5H7]. Parallel to 
these are statistical mechanics of DNA unzipping by force ]5]. In spite of a large amount 
of activity on the equilibrium statistical mechanics of these systems, there have been only 
few studies of the dynamics. The earliest study of the dynamics of such a system seems 
to be that of Sebastian [9] , who investigated the Rouse dynamics of a long chain molecule 
pulled out from a potential well. It was found that the escape of the polymer from the well 
happened if the pulling force / exceeded a critical value f c and that the time needed for 
escape of monomers scaled as n 2 d /{f — f c ). Similar results were obtained by Marenduzo 
et. al. who did simulations [10] . There have also been simulations of the desorption of a 
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pulled polymer [TT]. However, there has not been any analytic theory of the dynamics of 
this process, with excluded volume interactions included. In this letter, we develop such 
a theory based on the equations originally formulated by Brochard-Wyart et al [12] • The 
approach leads to a diffusion type of equation for the propagation of stretching along the 
polymer. As it accounts for excluded volume interactions, the equation is non-linear and 
in the mathematics literature it is referred to as diffusion equation with p-Laplacian |13j 
In the case of pulled dcsorption, one has to solve this non-linear equation with a moving 
boundary condition. Interestingly, an exact solution can be found for this problem for all 
the limits of interest. In the trumpet reigme, the time for desorption of n^ monomers r 
scales f c n 2 d / (j 2 ^ — fVj , very similar to what was found for the much simpler model 
of escape from a well by pulling [TJ]- If one adopted Rouse dynamics, this changes to 
t ~ fc^d/if ~ fc)- Surprisingly, in the other regimes too, time scales exactly in the same 
fashion as in the trumpet regime, though the dependence on the forces is changed. 

The Model. — We think of a polymer of N monomers, with excluded volume inter- 
actions. The unadsorbcd molecule would have a Flory radius Rp ~ aN v , where a is the 
size of a monomer. We shall label the monomer position along the countour length of the 
monomers with n which varies from — TV to 0. It is convenient to imagine this as a contin- 
uous variable. The solid surface on which the polymer is adsorbed is located at x = and 
at the time t = 0, one starts pulling on the end with n = with a force / in the positive 
x-direction. If the force on the monomer at the surface exceeds a critical value f c then the 
monomer is assumed to desorb. A simple view of the desorbed part of the molecule is in 
terms of the blob-picture introduced by Pincus [TSUI?] and used by Brochard-Wyart et. al 
[T2"1[TM2"0] and Sakaue [3T] to study the dynamics of pulled polymers. In this section, we 
shall discuss what has been referred to as the trumpet regime [3D], The other regimes are 
discussed in a later section. Let us denote the average position of the n'^monomer in the 
direction of pulling by x(n,t). As a result of pulling, the end at n = desorbs and starts 
moving as more and more monomers desorb. The desorbed part of the polymer may be 
imagined to consists of blobs of (varying) size R b streched out in the x— direction (see Fig. 
[1]). Denoting by An;, the number of blobs of size i?&,we can write 

x(n,t) = J2R b An b (1) 

6 

where the sum in Eq. |T]) is over all the blobs to the left of the n th monomer. Let us say that 
Ant blobs contain An monomers in total. As a blob of size Rb contains (Rb/a) 1 ^ monomers, 
where v is the Flory exponent, An and An;, are related through An;, = An/(Rb/a) 1 ^. 
Hence x(n, t) = a x l v J2b 1 ^An. The size of a blob is related to the force by Rb = ksT/ f, 

fcaT ) ' wnere 

—rid(t) labels the last monomer that is desorbed at the time t. On differentiating this with 
respect to n, we get 

This is the equation that we shall use in our analysis and is contained in the paper by 
Brochard-Wyart et. al [T2]. Using this we wish to calculate the average time of dcsorption 
of n th monomer. We now calculate the frictional force acting on the portion of the chain in 
the interval (x, x + Ax). Denoting the viscosity of the medium by r], a blob of size Rb would 
experience a frictional force rjvRb where v is the velocity with which the blob is moving. 
Hence the total force experienced by that region would be rjv ^2 b e & x Rb = ??|f Ax. Hence 
the net force acting on the n th monomer is / — r\ J dx ^ ^a"'^ dn. Using this in Eq. ©, 
we get the equation that describes the building up of tension along the chain on the right 
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hand side of the wall f!2l. 



dx 
dn 



k B T 



I 



dx dx 



(3) 



In Eq. ([3]), —Tid(t) < n < where n denotes any particular monomer on the desorbed part 
of the polymer. It may be noted that in the Rouse model, the analogue of Eq. ([2J would 
be linear in the force and may be obtained by putting v = 1/2. Further, each monomer 
experiences a force rjav (free-draining regime). Hence in the Rouse case, instead of Eq. 
we have 



dx 
On 



k B T 



I 



dx(n,t) 
Tja — — — dn 



(4) 




Fig. 1: The force / is applied in the positive direction, to the monomer at the end with n — 0, 
trying to desorb it. Desorption occurs if the force at the desorbing monomer (at n = —rid(t)) 
exceeds f c . x(n,t) denotes the average position of the n th monomer at the time t. The figure shows 
the stem-flower regime. In the trumpet regime, the stem is absent (i.e. n s (t) = 0). 



The Diffusion Equation with p-Laplacian . — It is convenient to use dimensionless 
variables (a tilde on top of the variable indicates that it is dimensionless) x = n = n/g, 
t = t/n with £ = k B T/f, g = (C/a) 1 ^, and n = r]^/k B T. This leads to 
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dx dx 
dt dn 



dn 



(5) 



Differentiating Eq. flSJ) gives ^ L 
becomes 



< 2 x 
'In- 



(II) 1_ " %• We sha11 P ut v = 3 / 5 and the equation 



d 2 x 2 fdx\* dx 



dt 



Putting n = in Eq. ((SJ) leads to 



dn 2 3 V dh 



dx 
dn 



(6) 



(7) 



n=0 



Further, at the monomer that is in the process of desorbing, the force has to have the critical 
value, which we denote by f c . Therefore, at the desorbing end, the boundary condition would 
be 

/i 1 -")/". (8) 



dn 



n=-n d (t) 
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In terms of dimensioned variables, f c = f c / f ■ That is, it is the ratio of the critical force to 
the actual force exerted, and in order to have desorption, f c < 1. The initial condition has 
all the monomers adsorbed. Thus 

x(n,0) = 0. (9) 

The equation (O is a diffusion equation with a p— Laplacian, with p = 3/2, which means 
fast diffusion in regions where ^ is small - this means that the tension spreads rapidly in 
regions where there is little tension 0. 

In the case of Rouse dynamics, Eq. ^ would lead to the diffusion equation ^| = ^0- 1^ . 
Thus with, or without the excluded volume interaction, the stretching of the polymer obeys 
a diffusion equation, and not a wave equation for which the disturbances propagate as pulses. 
Therefore, any movement of tension along the chain is diffusive. It is interesting to note 
that excluded volume interaction speeds up the diffusive motion (see ref [19]). 

The Solution. — As time passes, the polymer desorbs. We shall look for solutions of 
the form x(h,t) = t J h(u), u = and hd(t) = Vdt a . We take Vd to be a constant, to 

be determined. We need to choose values for the unknowns a, /3, 7 and v&. The boundary 
condition at the free end now reads 



dx 
On 



= p-f>h'(v d ) = l, (10) 

h=0 



which can be satisfied by putting 7 = /3 and imposing the condition 

h'{v d ) = 1. (11) 

Putting the form x(n,t) = t^h(u) into Eq. ([6]) leads to an equation for h(u) that has only 
the variable u in it if we choose a = /3 = 1/2. The resultant equation is 



3h"(u) + y/h'{u) ((u - Vd)h'(u) - h(u)) = 0. (12) 
In addition, the conditions © and © give 

h'(0) = ! 2 J\ (13) 

and 

h(0) = 0. (14) 

Eq. (|12p has to be solved subject to the initial conditions of (fT3"|) and Eq. ([TJJ. Further, the 
value of Vd is uniquely determined by the condition of Eq. ([lip . The numerically obtained 
x(n,t) is plotted against —n for different values of t in Fig. [5^. One can calculate the 
position of the pulled end of the polymer as a function of time, to get x(0,t) = sfihivd). 
In Fig. [5p we give a plot of Vd against the dimcnsionless force f c . We now analyze the 
behavior of Vd in the region where the pulling force exceeds the critical force only by a small 
amount. That is, f c ~ 1. In this case i>d <C 1 and one can do a perturbation analysis of Eq. 
(fT2"|) . We put h(u) = ho(u) + Vdh\(u) + Vjh2(u). Our aim is to evaluate h'(vd) upto terms 
quadratic in Vd- This can be achieved by evaluating hi(u) correct upto u 3_ \ Doing this we 
find that ho(u) = fc u, hi(u) = f c u 2 /6 and hz{u) = 0. Using these to evaluate h'{vd) gives 
h'{v d ) = fc/ 5 + fci>j/3. This can be used to solve Eq. ((TTJ) to get 



V3(l-/ c 2/3 )// c . (15) 
Thus the n, number of monomers that have desorbed at the time f is given by 

fid = v d f 1/2 - (16) 



1 An equation of the form §f = f^(^) p 2 ^ ls known as a diffusion equation with p-Laplacian. This 
has a diffusion coefficient proportional to 2 ( see reference |13l V 
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Fig. 2: (a) Plot of x(n,i) against — h for different values of t. (b) Vd against f c . The full curve is 
the result of numerical calculations while the dashed curve is the result of using Eq. (fT5)l . 



Defining r = ^ /kB T)vi-tafJk B T)^)) with T * = Va 3 /k B T, this becomes r = r n 2 in 
terms of dimensioned variables, for / close to f c . The case of Rouse dynamics was discussed 
in [14] . The result is (written in the notations of this paper) 

r = f c n 2 J (2C(/ - /„)) (17) 

for / close to / c , where £ = k B T jr\a? . This is very similar to the result with excluded volume 
interaction, except that the dependence on the force is different in the two cases. The result 
t ~ n 2 is in agreement with results of simulation, which found the total time of desorption 
of the molecule to be proportional to N 2 [11]. The above results show that as the force 
increases, excluded volume decelerates the desorption. This is a result of the non-linear 
stretching of the molecule with excluded volume, leading to an effective friction dependent 
on the pulling force. 

Larger forces: the stem and the stem-flower regimes . 

The Stem Regime ( f > f c > k B T/a). In this case, the molecule is stretched to the 
maximum extend, and would form what has been referred to as the stem |18j . The desorbed 
part has n s (t) monomers and a length of an s (t). As a result, it obeys r]a 2 n s (t)h s (t) = f — f c - 
Introducing the dimensionlcss variables n s (t) = n s (t), f = {af /k B T), f c = af c /k B T, and 
t = t(k B T/r]a 3 ), this may be solved to get n s (t) = \/2( f — f c )Vt, showing that in this case 
too, the desorption time scales like n 2 . 

The Stem Flower Regime ( f c < k B T/a < f) [H]. (see FigG}. The length of the stem 
region at any time can be found using the results of previous section, if one remembers that 
at one end of the stem there is a pulling force / and the other end has k B T /a. Therefore, the 
stem consists of n s (t) = \/2(/ — l)Vt monomers. The dynamics of the flower part is easily 
analyzed using the p— Laplacian equation. In terms of x = x/a, it reads ^§ = | (ff ) 2 
and has to be solved subject to the conditions (If)- - m = 1, (f^)- - fT , = fl^ and 

■> \ on I n=—n s (t) ' V cm / n= — n d (t) •"" 

x(n, 0) = 0. As earlier, putting x(n, t) = Vth(u), u = {n+nd(t)/Vt, and nd(i) = VdVt, leads 
to exactly the Eq. (fT2j). but with the conditions h(0) = 0, h'(0) = ft^ and h'(y c i — a) = 1 
where a = y/2(f — 1). The procedure of solution is just the same as earlier and the result 
is that t (x n 2 , the proportionality constant being known only numerically. 

Finally, we discuss the value of f c which depends on the adsorption regime. For the weak 
adsorption regime de Gennes et al. |22] have shown that the adsorption chain can be seen 
as a string of adsorption blobs with the size of each blob D oc a(e — e c )^ l/ ^ , where <f> is the 
crossover exponent. The dimensionless adsorption energy e = e/k B T and e c is the critical 
adsorption energy. The critical force results from the condition that the size of the Pincus 
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blob is equal to that of the adsorption blob, which leads to f c = ^T(e — e c ) v '^ /a. In the 
opposite limit of the strong adsorption we have f c = e/a > ksT/a. 

Summary and Conclusions. — We have analyzed the dynamics of desorption of an 
adsorbed polymer as a result of pulling at one of its ends with a force /. Our analysis takes 
into account the non-linear stretching of the polymer, which is due to the excluded volume 
interactions. The resultant p-Laplacian diffusion equation has to be solved subject to a 
moving boundary condition. Interestingly, exact solutions can be found in all the regimes 
of interest. In the trumpet regime we find that the time r required to desorb m monomers 
is given exactly by r = T nj with r = a 10 / 3 / c r// (3(fc B T) 4 / 3 (/ 2 / 3 - / c 2/3 )) , for / close to 

f c . For the simple Rouse model, this gets modified to r = f c n 2 d / (2£(/ — f c )). In the stem 
and stem-flower regime too, time has a quadratic dependence on the number of desorbed 
monomers, though force dependence is different. 

* * * 
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